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Abstract. Estimates of nuclear matrix elements for neutrinoless double beta decay (Oi^2/3) based 
on the quasiparticle random phase approximations (QRPA) are affected by theoretical uncertainties, 
which can be substantially reduced by fixing the unknown strength parameter gpp of the residual 
particle-particle interaction through one experimental constraint — most notably through the two- 
neutrino double beta decay (2i>2(3) lifetime. However, it has been noted that the gpp adjustment 
via 2i/2/3 data may bring QRPA models in disagreement with independent data on electron capture 
(EC) and single beta decay (/3~) lifetimes. Actually, in two nuclei of interest for 01/2(3 decay (^°°Mo 
and ^^^Cd), for which all such data are available, we show that the disagreement vanishes, provided 
that the axial vector coupling gA is treated as a free parameter, with allowance for < 1 ("strong 
quenching"). Three independent lifetime data (2^/2/3, EC, f3~) are then accurately reproduced by 
means of two free parameters {gpp, (7a), resulting in an overconstrained parameter space. In addition, 
the sign of the 2i^2/3 matrix element M^'^ is unambiguously selected (M^'^ > 0) by the combination 
of all data. We discuss quantitatively, in each of the two nuclei, these phenomenological constraints 
and their consequences for QRPA estimates of the 0v2/3 matrix elements and of their uncertainties. 
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1. Introduction 

The new paradigm of massive and mixed neutrinos, emerging from the recent evidence for neutrino 
flavor oscillations [H 121 [3l H], is still incomplete in several aspects. In particular, the nature of 
the neutrino fields (Dirac or Majorana) ^ remains undetermined, the amount of CP violation in 
the neutrino sector (if any) is unconstrained, and the absolute neutrino masses — as well as their 
ordering — are not yet known. The process of neutrinoless double beta decay (Oz/2/?), 

{Z,A)^{Z + 2,A)+2e- (0i/2/?) , (1) 

bears on all these issues and, thus, is a major research topic in current experimental and theoretical 
neutrino physics [6], [H El [9]. The claimed observation of Qv2j3 decay in ^^Ge with lifetime 
T^J^ ~ 2.2 X 10^^ y [in], and the projects aimed at its independent (dis) confirmation [9], have 
also given new impetus to the field. 

In general, barring contributions different from light Majorana neutrino exchange, the inverse 
Oz/2/5 lifetime in a given nucleus is the product of three factors, 

{T^,J,)-' = G"^ iM'^^l'ml,, (2) 

where G^" is a calculable phase space factor, M^'^ is the 0i/2/3 nuclear matrix element, and is the 
(nucleus- independent) "effective Majorana neutrino mass" which, in standard notation [11], reads 



1=1 



(3) 



iTii and Uf-i being the neutrino masses and the z/g mixing matrix elements, respectively. 

The calculation of the matrix element M^'^ in Eq. ([2]) for a candidate Oz/2/3 nucleus is notoriously 
difficult. It requires the detailed description of a second-order weak decay from a double-even 
"mother" nucleus (Z, A) to a double-even "daughter" nucleus {Z + 2, A) via virtual states (with any 
multipolarity J'^) of the so-called "intermediate" nucleus {Z + 1, A). The decay can proceed through 
both Fermi (F) and Gamow- Teller (GT) transitions, plus a small tensor (T) contribution, 

M"^' = M°"t + M"'^ - ^ , (4) 

9a 

and detailed nuclear models are required to estimate the separate components (X = F, GT, T). 
In the above expression, qa is the effective axial coupling in nuclear matter, not necessarily equal to 
its "bare" free-nucleon value qa — 1-25 |12j . 

Modern calculations of Oi^2/5 matrix elements are usually performed within either the 
quasiparticle random phase approximation (QRPA) [T3l [H] or the nuclear shell model (NSM) [15] 
and their variants, sometimes with large differences among the results. We remind that the QRPA 
basis of nuclear many-particle configurations, on which the residual particle-hole and particle-particle 
interaction is diagonalized to build the nuclear excitations, is limited to iterations of two-quasiparticle 
ones (reducing to the particle-hole configurations when the pairing interaction is switched off); for 
details, see, e.g. [H [7]. The advantage of the QRPA as compared to the NSM is that one can 
include essentially unlimited sets of single-particle states, even those forming the continuum of the 
positive-energy ones within the continuum-QRPA [16]. 
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Painstaking but steady progress in both the QRPA and the NSM approaches is gradually leading 
to a better understanding — and to a reduction — of the differences among their results [9]. However, 
even in the most refined approaches, the estimates of remain affected by various uncertainties, 
whose reduction is of paramount importance for both theory and experiment. Indeed, uncertainties 
in M^" propagate to the extracted value of (or hmit on) via Eq. ([2]), and affect directly the design 
of Qv2j3 experiments (in particular the detector size and the choice of the nucleus) needed to reach a 
given target sensitivity to m^/? ^j. Among the sources of uncertainties one can quote: (1) inherent 
approximations and simplifications of the theory; (2) existence of free or adjustable model parameters; 
(3) problematic description of the strong short-range repulsive interaction between nucleons; and (4) 
uncertainties in the value of qa- 

The latter problem arises from the significant reduction ( "quenching" ) of the strength observed 
in nuclear GT transitions (see, e.g., [l/Tj), which still lacks a clear experimental quantification and 
theoretical understanding. Two possible physical origins of the quenching have been discussed, one 
due to the A-isobar admixture in the nuclear wave function [18] and another one due to the shift of the 
Gamow- Teller strength to higher excitation energies induced by short range tensor correlations [T9] . 
In the absence of a better prescription, the effect of quenching (in either QRPA or NSM calculations) 
is often simply evaluated by replacing the bare value qa — 1-25 with an empirical, quenched value 
S'A — 1 [SO]- However, there is no a priori reason to exclude values (?a ^ 1, which have indeed 
sometimes been advocated, especially within the NSM approach p!5| [21]. 

In this context, we present a novel approach towards data-driven constraints on M^'^ calculations, 
assuming the possibility of strong quenching [qa < 1) within the QRPA. This unconventional 
hypothesis makes theory and data agree in a number of cases, where previous attempts have 
systematically failed. Therefore, we think that our approach may lead to a fruitful discussion and 
a fresh look at the whole problem of quenching, from both the theoretical and the experimental 
viewpoint. We stress, however, that we simply treat (7^ < 1 as a phenomenological possibility in this 
work, without any attempt to elaborate theoretical interpretations of the qa values emerging from 
the data analysis. 

Our work is structured as follows. In Sec. [2] we discuss the experimental data which can be used 
to benchmark the QRPA model. We adopt a selected data set, including the measured lifetimes of 
two-neutrino double beta decay, electron capture, and single beta decay for two nuclei, ^°°Mo and 
^^^Cd, which are of interest for searching 0^2(3 decay. In Sec. |3] we compare these data with the 
corresponding QRPA results, assuming standard quenching {qa = 1) or no quenching {qa = 1.25). 
We face then the well-known problem that the theory cannot match two or more data at the same 
time, for any given value of the so-called particle-particle strength parameter Qpp. In Sec. IHwe show 
that this problem can be phenomenologically removed if strong quenching [qa < 1) is allowed. In 
this case, the two parameters {gpp, qa) are overconstrained by three independent data in each of the 
two chosen nuclei, as shown in Sec. [51 In Sec. [6] we propagate the estimated [gpp, qa) uncertainties 
to the calculation of Oi/2/3 matrix elements and lifetimes, with and without the effects of short-range 
repulsive interactions. Finally, we summarize our results and discuss future perspectives in Sec. [3 



Overconstrained estimates of neutrinoless double beta decay within the QRPA 



4 



2. Experimental benchmarks 

In order to reduce the theoretical uncertainties, any nuclear model used in Oz/2/3 calculations should 
be benchmarked by as many weak-interaction data [22] as possible. Relevant weak processes are 
hsted in Eqs. ©-([Il]) below. 

Two-neutrino double beta decay (2z/2/?), 

(Z, A) ^ (Z + 2, A) + 2e- + 2Ve {2u2(3) , (5) 

is a second-order weak process (|AZ| = 2) which probes the same mother and daughter nuclei 
as Qv2(3 decay. It has been observed in several nuclei, thus providing a particularly important 
benchmark. Indeed, it was extensively demonstrated in |13] that the spread of QRPA calculations 
can be significantly reduced by constraining the nuclear model with the corresponding experimental 
21/2(5 decay lifetime (see [23] for earlier attempts). The 2v2(3 data help to fix an important free model 
parameter, namely, the strength Qpp of the residual particle-particle interaction [211 122] , and thus to 
"calibrate" the QRPA estimates of M^'^ . Despite the fact that the 2v2(3 decay process probes only 
a subset of the intermediate states relevant for Oz/2/3 decay (i.e., only those with J'^ = 1"^, via GT 
transitions), it is just the l"*" contribution to the total Oz/2/3 matrix element that reveals a pronounced 
sensitivity to Qpp, in contrast to the other multipole contributions [26]. This observation justifies the 
aforementioned fitting procedure employed in |13j . 

First-order weak processes (|AZ| = 1) related to 0i/2/? decay can probe, in usual jargon, either 
the "first leg" of the decay (from the mother nucleus to the intermediate one) or its "second leg" 
(from the intermediate nucleus to the daughter one). Relevant examples for the first leg include the 
electron capture (EC) from a bound state (e^J"), 

e,- + (Z + l,A)-^(Z,A) + z/e (EC), (6) 
and the charge-exchange reaction via (^He,t), 

3He+(Z,A)^(Z + l,A)+3H (3He,t), (7) 
as well as via (j9, n), 

p+{Z,A) ^ {Z + l,A) + n {p,n) . (8) 
The second leg is instead probed by the P~ decay, 

{Z + 1,A) ^ {Z + 2,A)+e' (/?-), (9) 
by the charge-exchange {d, ^He) reaction, 

2H + (Z + 2,A) ^ (Z + l,A) + 2He (d,2He), (10) 
and by ordinary muon capture (/iC), 

^- + {Z + 2,A)^iZ + l,A) + u^ ifiC) . (11) 

See also [27] for a recent discussion of these and other possible weak processes, including future 
(anti)neutrino-nucleus charged-current reactions at low energy [28]. Clearly, any of the above first- 
order weak processes could be used to set useful constraints on the nuclear model. Indeed, using 
P~ decay has been advocated as an alternative to 2z/2/3 decay for fixing the gpp parameter in QRPA 



Overconstrained estimates of neutrinoless double beta decay within the QRPA 



5 



Table 1. Compilation of experimental references for nine nuclear systems {A) of interest in Qv2[3 
decay ("^", "Z + 2" denote "mother", "intermediate", "daughter" nuclei, respectively). The 

entries refer to 2v2(3 decay data (|AZ| = 2) as well as to processes probing the so-called first and 
second leg (|AZ| — 1). For /9~ data, only decays from J'^ = 1+ states are considered. For muon 
capture (/iC), the data in 22] actually refer to natural isotopic mixture of the Z + 2 nucleus. See also: 
[111 115 for proposed EC measurements at A = 76, 82, 100, 116, and 128; |[44| for proposed (^He,t) 
measurements at A = 76 and preliminary (c?, ^He) data at ^ = 76 and 96; ^45] for preliminary /j,C 
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[29] [30| [3T] : /iC data might be similarly used in the near future [32]. However, one should be 
aware that these data are currently more sparse than for 2v2(3 decay and, sometimes, have inherent 
problems or limitations, as discussed below. 

Table [T] shows the current experimental status of the seven processes listed in Eqs. (ISHTTl). for 
nine nuclei of interest for Oz/2/5 decay searches. Data on 2z/2/5 decay lifetimes exist for all these nuclei 
[33] • Lifetimes for EC and (3~ decay have been measured only in three cases, A = 100, 116 and 128 
(with J'" = 1+ states for the intermediate nucleus) [311 [351 [361 [3Z]. In one case {A = 100), the most 
recent EC datum [35] appears to be in conflict with the older one [31] • Data on the charge-exchange 
scattering processes are also sparse. Available /iC data [12] are not particularly constraining, since 
they refer to the natural isotopic mixture containing the daughter nucleus; see however [16] for a 
comparison of QRPA calculations with /xC data, and [15] for preliminary /iC data in unmixed A = 76, 
82, and 150 daughter nuclei. Charge-exchange reactions involve analyses of spectral data which are, 
in general, more difficult to be interpreted and modeled than decay lifetimes |47l |48]. Data for 
(^He, t) exchange are available only for A = 100 and 116 [38]. In the latter case, the measured GT 
strength is in conflict with the one derived from EC [36j. Data for [d, ^He) exchange and A = 116 are 
reported in [H], where the GT strength distribution is, however, normalized to the reference (3~ one 
[37] at small excitation energy, and thus it does not provide an entirely independent constraint. The 
{p,n) reaction has been instead studied in several nuclei [39lllQ], with emphasis on the GT strength 
distribution (rather than on its normalization). For A = 116, it should be noted that the recent 
{p,n) data in [ID] disagree with the (^He,t) data in |38], and are only in rough agreement with the 
EC data in [36] . 

Clearly, new and dedicated measurements are needed, both to solve the mentioned experimental 
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Table 2. Experimental input. Half-life data (with la experimental errors) for 21/2(3, EC, and /3 
decay in ^"^Mo and ^^^Cd. All logarithms are in base 10. 



Nucleus 


log{T^J,/y) ± ae,p 
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discrepancies and to fill the missing entries in Tab. [T] [27l HH HI]. In the meantime, one needs to 
select a (hopefully consistent) data set, in order to perform a meaningful comparison with theoretical 
calculations. 

In this work we adopt the following approach: we ignore current data from the charge-exchange 
scattering processes (which, in several cases, either disagree with each other, or have no independent 
normalization, or provide poor constraints for our purposes), and we choose only those data which 
involve half-life measurements (rather than complex spectral analyses), namely, 2i>2(3, EC, and /3~ 
decay. Our investigation is then restricted to two nuclear systems for which all such data exist, 
namely, A = 100 and 116, which we shall often denote by the name of the "mother" nucleus (^°°Mo 
and ^^^Cd, respectively). For A = 100, we discard the old EC datum, log/t(EC) ~ 4.45 [53], in 
favor of the new (albeit unpublished) one [35] . Table [2] shows the corresponding input data that will 
be used in our analysis, in terms of logT/y (for Oz/2/3) and of log ft (for EC and where / is the 
usual nuclear Fermi function. (Throughout this paper, log = logj^g-) 

Although the (2z/2/5, EC, /3~) data are available also for A=128 (see Table 1), this nuclear 
system is left out of the consideration in the present work since the final nucleus ^^'^Xe is rather 
strongly deformed. The change in the deformation from an almost spherical ^^^Te to a rather well 
deformed ^^^Xe (/5 = 0.18 [49j) cannot be reliably treated within the spherical QRPA employed here. 
Importance of such an effect has been demonstrated in Refs. [501 EI] for the case of the 2z/2/3 decay 
using the deformed QRPA with schematic forces. 

3. Data versus theory with standard or no quenching (1 < gA < 1.25) 

In the context of the QRPA, it has been convincingly shown in [13] that the spread of theoretical 
calculations can be significantly reduced, in each of the nine nuclei in Tab. [H by fixing Qpp in such a 
way as to reproduce the measured 2Z/2/3 lifetimesfl This approach has however been questioned in 
[29t [30] . since the fitted value of gpp appears to underestimate (overestimate) the EC (/3~) lifetime 
by a large factor, as compared with experimental data. The alternative choice of fitting gpp by 
reproducing, e.g., the l3~ decay lifetime [301 [31], merely shifts the problem to other data (e.g., to the 
21/2/3 or EC lifetimes) which are no longer correctly reproduced; see also |52] for early examples of 
such a conflict. It is worth noticing that, in the related literature [TSj [31l [H], qa has been taken in 
the range 1 < (^a ^ 1.25, i.e., between standard quenching {qa — 1) and no quenching {qa — 1.25). 
Within such range, the problem of fitting two or more data (among 2v2f3, EC, I3~) appears to be 

X It is worth noticing that, in general, the effective values of both gpp and qa may change in different nuclei. 
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basically unsolved in the QRPA. Before discussing in the detail this problem in Sees. 13. 1H3. 21 below, 
we recall a few essential features of the QRPA. 

The (2i/2/3, EC, (3~) processes occur through GT transitions, either at first order in qa (for 
EC, I3~) or at second order in (for 2z/2/3). Therefore, theoretical estimates of the associated 
(logarithmic) lifetimes need to be performed only for = 1, and can then be scaled for ^ as: 



Within this work, QRPA calculations of the above lifetimes have been performed both in large 
basis (l.b., default choice) and in small basis (s.b.). The small basis consists of 13 single-particle 
levels (oscillator shells = 3 and 4, plus the f + h orbits from A^ = 5), while the large basis contains 
21 levels (all states from shells A^ = 1, . . . , 5), in accordance with the choice of [53l [T3] . The small 
set corresponds to 1^ particle- hole excitations, and the large one to about Ahuj excitations. 

An important output of QRPA calculations is the 2u2f3 matrix element M"^^ , whose modulus is 
probed by the observable Tl% according to 



where G^*^ is a calculable phase space factor, and the bare value of qa (1-25) is explicitly factorized 
out to make contact with previous notation [13]. In QRPA calculations, M^*^ typically starts positive 
for Qpp <^ 1, then decreases and eventually changes sign as Qpp increases. The critical value g*p 
where M'^^ = marks an infinite lifetime, logTy^ ~^ 00. It turns out that log ft{EC) is continuous 
across g*p, while log ft{j3~) diverges locally. For gpp increasing slightly beyond this critical point, 
the calculated energy of the first excited state Ei decreases and eventually vanishes, inducing a 
breakdown (the so-called "collapse") of the QRPA solution. QRPA calculations become thus less 
reliable in the vicinity of the critical and collapse points. 

Figure [U shows the matrix element M'^" as a function of Qpp for each of the two reference nuclei, in 
large basis. Similar results are found for small basis (not shown). In each panel, a vertical dotted line 
marks the critical value g*p where M'^'^ flips its sign. The value of M"^" drops rapidly for gpp > gpp, 
and the QRPA collapse is eventually reached. Both positive and negative values of M'^'^ may be 
phenomenologically acceptable in principle, although theoretical arguments suggest that M^'' > 
|13] . Determining the sign of M'^'^ is thus a relevant check of the theory. 

The QRPA estimates of M'^'^ , as well as those of the 2z/2/3, EC, and (3~ lifetimes, are affected 
by various sources of uncertainties. In Sec. Owe shall deal with the uncertainties related to the (a 
priori unknown) values of gpp and gA, and to the size of the basis. However, even if gpp and gA 
were perfectly known and the basis size were irrelevant, the approximation inherent to the QRPA 
approach would introduce further theoretical errors on each estimated lifetime. The assessment of 
these errors is obviously difficult and, to some extent, even arbitrary — but it is necessary to gauge the 
(dis) agreement between theoretical estimates and data. Our educated guess for the extra theoretical 
uncertainties (besides those related to gpp, to gA, and to the basis size) is ~ 20% for both the EC 
and 2z/2/5 lifetimes, and ~ 40% for the /3~ lifetime. In the latter case, a larger relative error is 
assumed, due to the smaller (by a factor 2-3) calculated values of the corresponding matrix element 



log(/t) ^ logif t/gl) for EC and /?" , 
log(Ti%) - log(Ti%/^7l) for 2Z.2/5 . 



(12) 
(13) 




(14) 
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as compared with the ones for the EC. Accordingly, we attach the following (±lo") theoretical errors 
o"th to each logarithmic lifetime, for any fixed values of ((^pp, qa) in any basis: 



In the next two subsections we shall compare the data in Tab. [2] with the corresponding QRPA 
estimates for = 1. It will be shown that, in none of the two reference nuclei, the QRPA results can 
be really made consistent with more than one datum at a time, within the quoted experimental and 
theoretical uncertainties. Moreover, it will become evident that higher values (e.g., gA = 1.25) 
can only worsen the situation. 

3.1. ^^^Mo data versus QRPA (gA = l) 

Figure [2] illustrates the comparison between ^^''Mo data and theoretical predictions for standard 
quenching {gA = 1) in large basis, as a function of gpp. The upper, middle, and lower panels refer 
to the 21/2/3, EC, and l3~ logarithmic lifetimes, respectively. In each panel, the horizontal band 
represents the experimental datum at ±lcr (as taken from Tab. [2]), while the curved band represents 
the QRPA results, with ilcr theoretical spread as in Eqs. ( |T5l - [T71) . Vertical dotted lines mark the 
critical value g*p which separate the left, positive branch (M^'^ > 0) from the right, negative branch 
[M'^" < 0). The preferred gpp ranges — where the experimental and theoretical bands cross each 
other — appear to be quite different in the three panels of Fig. [21 In particular, there is no overlap 
between the preferred gpp ranges in the upper and middle (or lower) panel, while there is only a 
marginal overlap between those in the middle and lower panels. Agreement between data and theory 
is never reached for all the three observables at the same time. 

If one choses the 2z/2/5 lifetime to fix ^fpp (as advocated in [13]), then two preferred ranges are 
selected, one in the positive branch (around gpp ~ 0.78), and the other in the negative branch (around 
^fpp ~ 0.79); see the upper panel of Fig. [2l Although both ranges are phenomenologically viable, the 
one in the positive branch is usually adopted on theoretical grounds [13] . However, for gpp — 0.78, 
the theoretical EC (/5~) lifetime turns out to be significantly smaller (larger) than the experimental 
value. Similar problems occur for gpp ~ 0.79 in the negative branch. 

Alternatively, one might use the /3~ lifetime to fix gpp (as advocated in [301 [3T ] ). In this case, as 
evident from Fig. [21 one could get marginal agreement between both f3~ and EC observables around 
gpp ^ 0.75, but only at the price of underestimating the measured 2i>2(3 lifetime by a factor of ~4. 
With one choice or another, it seems that current QRPA calculations fail to reproduce all the three 
independent lifetimes at the same time. 

The above discrepancies would become stronger by increasing the GT strength from its standard 
quenched value {gA — 1) to its bare value {gA — 1.25, not shown). For gA = 1.25, according to 
Eqs. ([T2l) and (fT3l) . the theoretical bands in Fig. [21 would be shifted downwards by —Hogg a — —0.4 
(upper panel) or by — 21og5'A — —0.2 (middle and lower panels). The preferred ranges of gpp would 
then move to the right for 2z/2/5 and f3~ , and to the left for EC, thus destroying even the marginal 
agreement existing between j3~ and EC observables for gA = l- We conclude that, within the range 



log(T,%/y) 
log/t(EC); 
log/t(r): 



att = ±0.08 , 
ath = ±0.08 , 
ath = ±0.15 . 



(15) 
(16) 
(17) 
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1 ^ 9a ^ 1.25, current QRPA calculations cannot reproduce the three hfetime data (nor, to some 
extent, any two among them) for any value of Qpp. These graphical results will be numerically 
confirmed in Sec. [5l 

3.2. ^^^Cd data versus QRPA (gA = 1) 

Figure [3] is analogous to Fig. [2], but for ^^^Cd. The situation is very similar to ^"'^Mo, and the same 
qualitative considerations apply, although the preferred ranges of gpp are different. Also in this case, 
it is not possible to reconcile the QRPA estimates with the three independent lifetime data data for 
any value of Qpp, at fixed qa = ^- The discrepancy becomes worse for qa = 1.25 (not shown). 

We remark that Figs. [2H3] refer to QRPA calculations in large basis. Very similar results are 
obtained — and the same comments apply — to calculations in small basis (not shown). 

4. Data versus theory with strong quenching (gA < 1) 

In the previous Section, we have shown that the QRPA fails to reproduce the three (2i/2/3, EC, /?") 
lifetimes in each of the two reference nuclei (^°°Mo and ^^^Cd), as far as qa is taken in the usual 
range, 1 < < 1.25. In particular, the discrepancy becomes worse as one moves towards the upper 
end of this range. Conversely, the discrepancy can be expected to become less severe (and hopefully 
vanish) for qa < I, corresponding to a "strong quenching" of the GT coupling. 

Values of qa lower than unity, although rather unconventional in the QRPA literature, are 
not uncommon in NSM calculations. The NSM, being an ab initio approach, does not depend on 
phenomenological parameters such as Qpp, but of course retains the dependence on the axial coupling 
gA, with the associated quenching uncertainties. Although a quenched value (/^ ~ 1 seems to 
roughly provide the correct normalization of the GT strength, strongly quenched values qa < ^ may 
occasionally be needed to bring NSM calculations in agreement with data [151 HI] . It is fair to say 
that, in the NSM approach, one is not committed to a strict range for qa (such as 1 < (^yi < 1.25): 
any value qa ~ 0(1) is generally accepted, if the data require so. 

In both the QRPA and the NSM approach, the origin and size of the GT quenching remains 
in part obscure and uncertain from a theoretical viewpoint, and the inferred values of qa fluctuate 
considerably in different data analyses, processes, and nuclei. Even for a fixed process and nucleus, it 
is not excluded that the quenching may be energy- dependent [46]. Therefore, the common practice of 
adopting either the standard quenched value qa — ^ or the bare value qa — 1.25 may be unnecessarily 
restrictive. It is perhaps more sensible to treat qa as a free parameter of order unity, whose precise 
value needs to be constrained by the data themselves, rather than pre-assigned by theory — ^just as 
one does for gpp. In the following, we thus adopt a purely phenomenological viewpoint, and show 
that specific choices of qa below unity (which will be more precisely derived in Sec. [5]) can bring 
QRPA calculations in agreement with all the three lifetime data in each of two reference nuclei. 

4.1. ^°^Mo data versus QRPA (gA = 0.74) 

Figure H] is analogous to Fig. [21 but for the strongly quenched value qa = 0.74. We anticipate 
that this value provides the best overall agreement of QRPA calculations (curved bands) with the 
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data (horizontal bands). Around Qpp ~ 0.73, all bands cross each other in the three panels. No such 
common crossing occurs in the negative branch, as also confirmed by numerical explorations. Besides 
selecting the positive branch, the data appear to prefer a particle-particle strength {gpp ~ 0.73) 
sufficiently far from the the critical and collapse values, where the QRPA estimates become less 
reliable. The [3~ theoretical band in the lower panel is rather steep around Qpp ~ 0.73, and can thus 
provide, together with the experimental datum, both an upper and a lower bound to Qpp, the upper 
(lower) bound can also be enforced by the 2v2(3 (EC) observable, as evident in the upper (middle) 
panel. In perspective, a reduction of the EC error in the middle panel would be beneficial to better 
probe this strong-quenching scenario. 

4.2. ^^^Cd data versus QRPA (qa = 0.84) 

Figure [5] is analogous to Fig. [3l but for the strongly quenched value qa = 0.84. A good overall 
agreement between theory and data is reached in a broad range Qpp ~ 0.4-0.6. It is interesting to 
note that this range could be significantly restricted if the experimental errors of the EC datum [36] in 
the middle panel were reduced by a factor of two or more. Also in this case, the data unambiguously 
select the positive branch, and keep Qpp far from the critical and collapse points. 

We remark that Figs. HHS] refer to QRPA calculations in large basis. Very similar results are 
obtained — and the same comments apply — to calculations in small basis (not shown). 

4.3. Discussion 

Strong quenching {qa < 1) appears to provide a phenomenological solution to the well-known overall 
discrepancy between QRPA results and lifetime data. This solution is nontrivial because: (1) two 
free parameters enable to reproduce very well, within la uncertainties, three independent data (in 
each of two different nuclei); (2) the positive branch {M'^" > 0), which is favored by theoretical 
arguments, is unambiguously selected by the data; (3) the preferred values of gpp are far enough from 
the critical and collapse values. Such data-driven features seem to be more than accidental facts, 
and suggest that qa < 1 might be a realistic option within the QRPA. More accurate lifetime data 
(especially for EC and, to some extent, for 21/2/3 decay), as well as further charge-exchange reaction 
data (not considered in this work) should provide additional probes of the strong quenching solution. 

This solution is admittedly unconventional in the context of QRPA, where qa has been 
customarily taken within the range 1 ^ Qa ^ 1.25. It may be that strong quenching is associated to 
other effects, whose degrees of freedom might be traded for milder variations of qa- However, if new 
free parameters are added to gpp and qa, the data set must also be enlarged to provide meaningful 
and nontrivial constraints — not much would be learned, in general, by fitting data with > 
parameters. 

For the sake of simplicity, in this work we do not explore more elaborate scenarios with additional 
data and further QRPA degrees of freedom. We just take for granted the indication in favor of 
qa < 1, and perform quantitative fits to three selected data (the 2i/2/?, EC, and P~ lifetimes) via two 
parameters {gpp, qa)- We shall thus obtain an overconstrained parameter space, used for subsequent 
Oz/2/3 calculations in Sec. [61 Despite the above caveats, this approach represents a step forward with 
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respect to previous attempts, which aimed at reducing the Oi/2/3 model uncertainties in QRPA by- 
fitting a single datum (either 2z/2/5 or (3~) through a single parameter {gpp) at fixed qa- 

5. Overconstraining the {gpp, qa) parameters 

We perform a least-square fit to the three data xi = \og(Ty2/j)j ^2 = log /t(EC), and 0:3 = log ft{P~) 
in terms of the two free parameters {gpp, Qa) ■ The function to be minimized is defined as 



where all the ingredients have been defined in the previous Sections. Asymmetric experimental errors 
(see Tab. |2]) are properly included by choosing either the upper or lower error, according to the sign of 
the difference xf^ — xf^'^. The minimum search is performed by numerical scan over a dense grid in the 
{gpp, Qa) rectangle [0, 1] ® [0, 1.25]. Given three data and two parameters, one expects Xmm ~ 
for a proper fit. The expansion around the best-fit values of {gpp, qa) at Ax^ = ~ Xmin = "-^ 
provides then the n-a contours for such parameters [iTj. In the following, we show the main results 
both in graphical and tabular form. 

Figure [6] shows the results of the {gpp, gA) fit in large basis. In each of the two panels 
(corresponding, from top to bottom, to ^°°Mo and ^^^Cd) a dot marks the best-fit point, surrounded 
by the 1, 2 and Sa contours. Vertical dotted fines separate the positive and negative branches of 
M^*^. In both panels, the allowed regions are fully contained in the positive branch, thus confirming 
quantitatively the theoretical arguments in favor of M'^'^ > [13]. The best-fit points are safely 
far from extremal values of gpp (0 and ^ gpp), but the allowed regions may extend towards one of 
them. In particular, the allowed range of gpp is somewhat squeezed towards the critical value for 
^°°Mo, while it extends towards zero for ^^^Cd at 3a. More accurate experimental data (especially 
from EC and, to some extent, from 2i^2/5 decay) would be useful to shrink such ranges, as discussed 
in Sec. HJ and might thus prevent the occurrence of nearly extremal values of gpp. Concerning gA, 
strong quenching {gA < 1) is definitely preferred at > 3a in both cases. 

We emphasize that "overconstraining the {gpp, gA) parameters" is equivalent to state that, 
in each of the ^°°Mo and ^^^Cd reference nuclei, our scenario makes one prediction which is 
experimentally verified. Figures [7] and [8] illustrate this statement via the la bands individually 
allowed by (3~ , EC and 2v2f3 data for ^''^Mo and ^^^Cd, as obtained by a breakdown of the three 
contributions in Eq. (fT8|) . Any two bands can be used to constrain {gpp, gA) in a closed region (the 
"prediction"), which is then crossed by the third independent band (the "experimental verification"). 

The numerical results of the global {gpp, gA) fit in large basis are summarized in Table [3l The 
fit quality is very good in all cases (Xmm ~ 1) ^^"^ best-fit values for the three lifetimes are in 
striking agreement with the corresponding data in Tab. [21 which are repeated for convenience in 
Table |3] (in square brackets). The best-fit values and incr ranges (n = 1, 2, 3) for gpp and gA are 
also reported. (The gA values adopted in Figs. [Hand [5] are just taken from Table [3l) 

We have repeated the analysis in small basis, with similar results. The graphical results are 
omitted, while the numerical ones are reported in Table |H The quality of the fit is very good also 
in this case. The allowed ranges for gpp and gA in small basis (Tab. H]) are somewhat different from 




(18) 
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Table 3. Results of the (gpp, qa) fit for two different (mother) nuclei, with QRPA calculations 
performed in large basis. Column 2: minimum x^. Columns 3-5: theoretical lifetimes for 2v2j3, EC 
and j3~ decay at best fit, to be compared with the experimental data in Tab. [2] which are repeated 
here in square brackets. Columns 6-9: value of gpp at best fit, and allowed ranges at f , 2 and 3(t. 
Columns 10-13: value of at best fit, and allowed ranges at 1, 2 and 3cr. 



Nuclei 


Xmin 


\og{TljJy) 


log/i(EC) 


log/t(r) 


9pp 


±1(7 


±2(7 


±3(7 


9A 


±1(7 


±2(7 


±3(7 


looMo 

"6Cd 


1.26 
0.12 


18.82 [18.85] 
19.49 [19.48] 


4.09 [3.96] 
4.35 [4.39] 


4.66 [4.60] 
4.63 [4.66] 


0.733 
0.493 


+0.020 
-0.020 

+0.106 
-0.149 


+0.031 
-0.063 

+0.173 
-0.358 


+0.039 
-0.126 

+0.224 
-0.493 


0.741 

0.843 


+0.046 
-0.037 

+0.042 
-0.037 


+0.120 
-0.074 

+0.088 
-0.075 


+0.176 
-0.107 

+0.149 
-0.106 


Table 4. As in 


Table H but 


in small basis. 




















Nuclei 


■ 

Amin 


\og{TllJy) 


log/i(EC) 


log/t(r) 


9pp 


±1(7 


±2(7 


±3(7 


9A 


±1(7 


±2(7 


±3(7 


looMo 

"6Cd 


1.11 

0.03 


18.82 [18.85] 
19.49 [19.48] 


4.08 [3.96] 
4.37 [4.39] 


4.67 [4.60] 
4.65 [4.66] 


0.862 
0.540 


+0.024 
-0.035 

+0.130 
-0.165 


+0.043 
-0.094 

+0.220 
-0.385 


+0.055 
-0.181 

+0.283 
-0.538 


0.745 
0.815 


+0.042 
-0.037 

+0.042 
-0.033 


+0.098 
-0.074 

+0.084 
-0.070 


+0.172 
-0.111 

+0.139 
-0.102 



those in large basis (Table [3l), but with similar features. In particular, the allowed Qpp range is in the 
positive branch, and the general trend in favor of qa < ^ is confirmed. We conclude that the main 
results obtained so far do not change qualitatively with the size of the basis. 



6. Implications for Oz/2/3 decay 

In the previous Section we have obtained allowed regions in the parameter space [gpp, Qa)- In this 
Section we study how such regions affect the QRPA calculation of Oz/2/5 decay, after recalling some 
basic features of this process. 

The (2z/2/3, EC, processes that we have considered so far occur only via GT transitions 
through 1+ intermediate states. The leading contribution MqI^ to the amplitude of the neutrinoless 
double beta decay also comes from the GT-type transitions which, however, proceed through 
intermediate states of all, but 0^, multipolarities. In addition, there are Fermi M^" and (small) 
tensor M^y contributions to the Oi/2/3 matrix element, 

M'-'igpp, qa) = M'c^igpp) + M^^gpp) - ^ (19) 

9a 

where the dependence on gpp and gA is made explicit. 

Figure [9] shows the relevant components of the Oz/2/5 matrix elements as a function of gpp in 
large basis, and including short range correlations, which will be shortly discussed below. Since the 
QRPA calculation is computer-intensive, gpp is varied only within the relevant ±3o" range shown in 
Fig. [6l Note that the leading component shows significant variations with gpp, so that any constraint 
on this parameter (such as those derived in the previous Section) helps to reduce the spread of 
QRPA estimates of Oz/2/3 decay. Results qualitatively similar to Fig. [9] are obtained for small basis, 
or without short range correlations (not shown). 
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Given the QRPA results in Fig. [HI the Oi/2/3 matrix element can be computed for any relevant 
value of Qa and Qpp through Eq. f|T9|) . In order to make contact with the notation in Ref. [13j, we 
shall actually rescale the matrix element as 



M'O" = M"'^ (^^j . (20) 
The Oi/2/3 lifetime reads then 

T^''{9pp,9a) = I ^,0.^1 2 , (21) 
where the proportionality factor t'j'/2 = (m|^G'°'^/1.25^) ^ [y] is numerically given by 

,0. _ / 1-83 X 10^^ rUo) , 

- \ 1.68 X 1027 (ii^Cd) , ^ ' 



for a reference Majorana mass m^^ = 50 meV. For different values of m^^, one just rescales 

^1/2 ^I30- 

For any given value of [gpp, Qa), calculations of M'^" are affected not only by the size of the basis 
(either large or small), but also by uncertainties which are peculiar of the 0i/2/3 process, namely, 
those related to the important issue of short range correlations (s.r.c). These correlations account 
for the well-known fact that the nucleon-nucleon interaction becomes strongly repulsive at small 
internucleon distances. This in turn must lead to strong suppression of the relative-motion wave 
function at small distances (s.r.c. effects). Short range correlations are explicitly included neither 
within the QRPA nor within the NSM. They are instead introduced ad hoc directly into the neutrino 
potential via a multiplicative factor (the square of a correlation function). One of the most popular 
is the Jastrow-like correlation function [51] which has been used in the previous calculations [53l [13] 
and is also used in this work. We shall thus present results in four cases, corresponding to either 
large or small basis, with or without the Jastrow-like s.r.c. effects. 

In each of the four cases, the effect of the {gpp, qa) uncertainties on M'^^ is estimated by 
marginalization [IT] , taking into account the fact that the same fixed value for the matrix element may 
be realized by different ("degenerate") couples of values {gpp, qa)- More precisely, given the function 
X^{9ppy 9a) defined in the previous Section, and for a fixed value M'^" , we define a marginalized 
function, 

X'(M'°'^)= minx'(^PP,^A) , (23) 
over the degenerate set of {gpp, (Ja) obeying 

M"^{~gpp,~gA)=M'"' . (24) 

The minimization of x'^{M'^'^), and the expansion around the minimum at Ax^ = n^, provide the 
correct best-fit values and na ranges for M'^^ , respectively. Since we are interested in < 3, we 
perform a numerical marginalization over a dense, rectangular grid covering only the ±3cr ranges of 

{gpp-, gA)- 

Tables [5] and [6] provide an overview of the derived ranges for M'^^ at 1 , 2 and 3(T (in large and 
small basis), with and without the effect of s.r.c, respectively. We also report the corresponding 
ranges for the measurable (log) lifetime T^^^g, at the reference value mp = 50 meV. Note the ±na 
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Table 5. QRPA estimates of the Qv2(3 matrix elements M"^'^, including the effect of short range 
correlations. The central value and the allowed ranges of M"^'^ are derived, respectively, from the 
best-fit values and allowed ranges of {gpp,gA)- The estimates refer to both large basis (l.b.) and 
small basis (s.b.). We also report the corresponding (logarithmic) ranges for the Qi'2/3 lifetime rjy'2 
(in years), assuming a reference value mpp — 50 meV. 



Nucleus 




±1(7 


±2a 


±3(7 


iog(r{';2) 


±lcr 


±2cr 


±3cr 


Large basis 


















looMo 


2.66 


+0.15 
-0.14 


+0.33 
-0.25 


+0.61 
-0.35 


26.411 


+0.046 
-0.046 


+0.088 
-0.100 


+0.124 
-0.180 


"6Cd 


2.44 


+0.23 
-0.18 


+0.53 
-0.32 


+0.90 
-0.44 


26.448 


+0.065 
-0.079 


+0.123 
-0.169 


+0.174 
-0.272 


Small basis 


















looMo 


2.45 


+0.16 
-0.15 


+0.35 
-0.25 


+0.65 
-0.34 


26.485 


+0.055 
-0.055 


+0.095 
-0.116 


+0.132 
-0.204 


"6Cd 


2.15 


+0.20 
-0.16 


+0.46 
-0.29 


+0.78 
-0.41 


26.561 


+0.067 
-0.079 


+0.127 
-0.169 


+0.181 
-0.271 



ranges are generally asymmetric and do not scale linearly, in part as a consequence of the original 
one-sided gpp limits at either or ^ g*p (see Fig. E]). By comparing the results in Tables and El 
it appears that the basis size is not the major source of systematic uncertainties. Conversely, the 
inclusion or exclusion of s.r.c. effects always induce changes > la. 

Figure [10] shows an overview of QRPA results for the nuclear matrix elements (including s.r.c. 
effects) in three different cases for each nucleus. From left to right, the first two cases correspond to 
the la ranges from Table [5], in large and small basis, respectively. The third case correspond to the 
results previously obtained in [13] for qa = 1 (with correspondingly smaller error bars, due to the 
fixed Qa value). Remarkably, such results for M"^^ p[3] differ by < 12% from those obtained in this 
work, in spite of a marked difference in the central values of qa and gpp. 

Summarizing, in each of the two nuclei examined it is possible: [i) to fit very well three data 
(2i/2/5, EC, /3~) with two parameters {gpp,gA), provided that gA < 1; {ii) to exclude the negative 
branch M"^" < 0; and (Hi) to derive robust ranges for 01/2(3 observables. There remains a relative 
large uncertainty on the 0i>2(3 matrix element, associated with the size of short range correlation 
effects. Unfortunately, s.r.c. effects are peculiar of 0^2(3 decay and are not constrained at all by the 
(2i/2/5, EC, P~) data considered in this work. 

7. Conclusions and Perspectives 

It was shown in [13\ that, by fitting gpp in order to reproduce in calculations the corresponding 
experimental 2i>2f3 decay lifetimes, the sensitivity of calculated Oz/2/3 matrix elements to other 
ingredients of the QRPA, such as the basis size, can be successfully removed. Also, it was shown 
that the sensitivities of the results to gA gets much milder than one could naively expect. There are 
also different proposals for fixing gpp, for instance, by reproducing the single beta decay observables 
as advocated in [29l [30]. By fitting gpp to reproduce the P~ lifetimes of the ground states of the 
intermediate nucleus one gets the results which are similar to the ones obtained in fl3\ , but the EC or 
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Table 6. As in Tab.[5l but without short range correlations. 



Nucleus 




±lcr 


±2cr 


±3cr 




±lcr 


±2(7 


±3(7 


Large basis 


















looMo 


3.27 
2.84 


+0.16 
-0.15 

+0.25 
-0.19 


+0.34 
-0.29 

+0.57 
-0.35 


+0.63 
-0.42 

+0.98 
-0.49 


26.233 
26.317 


+0.041 
-0.043 

+0.061 
-0.073 


+0.080 
-0.087 

+0.116 
-0.159 


+0.117 
-0.155 

+0.165 
-0.256 


Small basis 


















looMo 

"6Cd 


2.97 
2.47 


+0.17 
-0.16 

+0.22 
-0.17 


+0.36 
-0.29 

+0.49 
-0.32 


+0.67 
-0.40 

+0.84 
-0.44 


26.318 
26.440 


+0.048 
-0.049 

+0.063 
-0.075 


+0.089 
-0.100 

+0.119 
-0.159 


+0.125 
-0.178 

+0.171 
-0.255 



2v2(3 lifetimes are not reproduced. In this paper we have tried to reconcile all these data (available 
for the two nuclei ^°°Mo and ^^^Cd) by letting qa to be a free parameter of the model. In each 
nucleus, we have then found systematic indications in favor of strong quenching [qa < 1), and we 
have been able to overconstrain two parameters ((^pp, qa) with three lifetime data (2z/2/3, EC, 
as well as to fix the sign of M"^^ (> 0). 

The quenched values of qa for A = 100 and A = 116 nuclear systems obtained in this work 
(qa — 0.74 and qa — 0.84, respectively), although a bit unusual, are not much below the typical 
range qa — 0.9-1.0 (corresponding to the quenching factor q = g a/ 1-25 ~ 0.7-0.8) used within the 
NSM for lighter nuclei [I5]. Even stronger quenching q ~ 0.5 (corresponding in our notation to 
qa ~ 0.6) has been called for in shell model calculations [55| |56| [57] of the Gamow- Teller strength 
for nuclei in the region of A ~ 100, to which the systems considered in the present work are close. 

The physical origin of the quenching of qa has been discussed in the past. One explanation [18] 
assigns this effect to the A-isobar admixture in the nuclear wave function. Another — more generally 
accepted — explanation [19| assigns the quenching to the shift of the Gamow- Teller strength to higher 
excitation energies due to the short range tensor correlations. In light nuclei the quenching found in 
Ml transitions reduces qa from its bare value (~1.25) to the in-medium one (~1). But the actual 
quenching in nuclear structure calculations can depend as on the detailed nuclear environment as 
on the truncations inherent to the model such as, for example, the basis size. Therefore, it appears 
useful to revisit the theoretical explanations of quenching, in order to check if and how they can 
cover cases with (7^ < 1, as those emerging from our phenomenological analysis. 

From the experimental viewpoint, it has been already mentioned that future EC data [l3l HI] 
will be especially relevant in improving the {gpp, Qa) parameter constraints. Moreover, the strong 
quenching of the axial vector coupling constant qa should be observed not only in single and double 
beta decays, but also in Ml transitions. Therefore, the study of charge-exchange reactions as (p, n), 
(n, p), (^He,t) and (rf, ^He) [2ll [41], |48] can shed new light on this issue. It is imperative, however, 
that the data are analyzed with no prior or hidden hypotheses about the GT coupling qa- 

In conclusion, we think that the results of this work offer a novel possibility to reconcile 
QRPA results with experimental data, which deserves further discussions and tests, and warrants a 
revisitation of the quenching problem from a new perspective. By the present analysis, we are able 
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to assign in a controlled manner theoretical uncertainties to the calculated matrix elements for the 
Oz/2/9 decay. Remarkably, our present results for M'^'^ agree within the error bars with those obtained 
in [13] for qa = 1.0. 
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Figure 1. Matrix elements for 2v2(3 decay in the QRPA (solid curves) as a function of g^p, for 
iooMo and ^^^Cd. In each panel, the vertical dotted line marks the critical Qpp value where M21, = 0. 
Calculations refer to the large basis. 
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Figure 2. Lifetimes for the 2z^2/3, EC, and /3~ decay in ^'^"Mo. Horizontal bands: experimental 
data. Curved bands: the QRPA results as a function of gpp, for fixed = 1, in large basis. The 
vertical width of the bands corresponds to ±la uncertainties. 




Figure 3. Lifetimes for the 2v2j3, EC, and /3~ decay in ^^^Cd. Horizontal bands: experimental 
data. Curved bands: the QRPA results as a function of gpp, for fixed = 1, in the large basis. The 
vertical width of the bands corresponds to ±la uncertainties. 
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Figure 4. As in Fig. [51 but for gA = 0.74. Note the overall agreement of the QRPA results with the 
data for gpp ~ 0.73. 
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Figure 5. As in Fig. [31 but for gA — 0.84. Note the overall agreement of the QRPA results with the 
data for gpp ~ 0.5. 
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Figure 6. Regions allowed at n-a in the {gpp, qa) plane from a QRPA fit to the 2i/2f3, EC, and /3 
data, in each of the two nuclei ^"^Mo and ^^®Cd. The QRPA calculations refer to the large basis. 
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Figure 7. Breakdown of individual constraints in the {gpp, gA) plane for ^"°Mo. The slanted bands 
corresponds to the regions allowed at la level (including experimental and theoretical errors) by /3~, 
EC, and 2v2[3 data. Their combination (thick ellipse) coincides with the \a contour in the upper plot 
of Fig. 




Figure 8. As in Fig. [71 but for ^^^Cd. The thick elhpse coincides with the la contour in the lower 
plot of Fig. M 
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Figure 9. Oi^2/3 matrix element components A/q^^ + M^'^ (solid) and —Mp"^ (dashed), as a function 
of the gpp parameter in its 3a allowed range (see Fig. [6]). The QRPA calculations refer to the default 
case (the large basis with the Jastrow-like short range correlations). 
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Figure 10. Overview of Oi^2/3 matrix elements Af'"'^, togetlier with their ztlcr estimated errors. For 
each nucleus, three QRPA cases are shown. From the left to right, the first two cases correspond to 
the results of this work in the large basis (black circle, with thick error bars) and in the small basis 
(black square). The error bars for these two cases encompass the uncertainties in both parameters 
{Qpp, gA) from the fit to {2i/2f], EC, /3~) data. The third case (white circle) refers to the previous 
results of Ref. [13], as obtained for the fixed value gA = ^ (with gpp adjusted to 2i/2/3 data). All cases 
include the effects of the Jastrow-like s.r.c. 



